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Two-dimensional QED with N flavor fermions is solved at zero and finite tem- 
perature with arbitrary fermion masses to explore QCD physics such as chiral 
condensate and string tension. The problem is reduced to solving a Schrodingcr 
equation for TV degrees of freedom with a specific potential determined by the 
ground state of the Schrodinger problem itself. 



1 QCD 4 physics in QED 2 

Two-dimensional QED with massive fermions is not exactly solvable, ie. it 
is not reduced to a free theory like the massless case. When the fermion 
masses are large it becomes a highly interacting model. Due to its many 
similarities with four-dimensional QCD in such respects as instantons, chiral 
dynamics, and confinement, it has long been a testing ground for intuitive 
physical principles and new ideas. 

One of the purposes of this work is to clarify these aspects of QCD physics 
by evaluating chiral condensates, the Polyakov loop, and string tension as 
fermion masses (to), vacuum angle (8), and temperature (T) vary. We shall 
show that N flavor QED is effectively reduced to the quantum mechanics of 
N degrees of freedom, which can be solved numerically on workstations Era 

The Lagrangian is 

N 

C = -lF fl „F» v + Y,Tpa{7»(id fl -eA l t)-m a }ij a . (1) 



a To appear in the Proceedings of Continuous Advances in QCD 96, University of Min- 
nesota, March 28-31, 1996. 



1 



We examine the model denned on a circle S 1 with a circumference L. Upon 
imposing periodic and anti-periodic boundary conditions on the bosonic and 
fermionic fields, respectively, the model is mathematically equivalent to a the- 
ory defined on a line (R 1 ) at finite temperature (T) by analytic continuation 
to imaginary time r (= it) and interchange of r and x. Various physical quan- 
tities at T ^ on are obtained from the corresponding ones at T = on 
S 1 by substituting T _1 for L. 

2 Zero-mode Hamiltonian 

To begin the computation we bosonize the fermions. Each two-component 
fermion (ip a ) is expressed in terms of zero modes {qtiPt) and oscillatory 
modes (<p a (x), n a (x)). The boundary conditions enforce that pf 's take integer 
eigenvalues. p~ — p+ and p~ + p+ correspond to charge and chiral charge, 
respectively. In a vector-like theory without background charges one can stay 
in a subspace defined by p^ = p~ . 

The relevant parts of the Hamiltonian are expressed in terms of {q a = Qa + 
<laiVa=\[vt +Pa}), (0a,IT Q ), and (8 W , Pw) where 6 W is the Wilson line 
phase, the only physical degree of freedom associated with gauge fields on a 
circle. The Hamiltonian becomes 

Htot — Ho + H^, + -ffmass 
It i N 

a=l 

u v a—1 v a ' 

where /i 2 = Ne 2 /it. H mass represents the contribution coming from the fermion 
mass terms. Expression (||) is exact, from which it immediately follows that 
QED2 with massless fermions is exactly solvable. It contains 1 massive boson 
and N — 1 massless bosons. 

The fermion masses give nontrivial interactions among the zero modes and 
4> modes, and the previously massless bosons now become massive. To find the 
true vacuum, we first determine the vacuum wave function in the zero mode 
sector with given physical boson masses /i a 's (a=l,- • • , N). The boson masses, 
which depend on the vacuum structure in the zero mode sector, are recomputed 
with the vacuum wave function thus obtained. Since input and output values 
for the boson masses must be the same, this gives a self-consistency condition 
which we can solve for numerically in general, and analytically in certain limits. 
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As a basis one may take eigenstates of (PwiPa) with eigenvalues (j>w, n a) 
where the n a 's are integers. It is more convenient however to take a coherent 
state basis with respect to {n a }. The vacuum wave function is written as 
fipwl fi, • • • i <Pn—i> One of the angular variables, 0, specifies the so-called 
^-vacuum; due to gauge invariance its value does not change. 

When fermion masses are small, / = e _7r ' iip »'/ 2Jv y((^; 9) to good accuracy. 
/(</?; 9) must satisfy 

| - A N + Vn\ f(fir ■ ■ , tpN-i) = e f(<pi, • ■ - ,<pn-i) 

s 2 2 ^ a 2 



A \ " J± j_ \ ' ^ 

w ^<v 2 N-\j£d<p a d<p h 

N N 

V N = - ^2 m a A a COS ip a (^2 i Pa=&j (3) 



where A a is determined by the boson masses ji a . The problem is now a 
Schrodingcr equation. The salient feature is that the potential has to be de- 
termined self-consistently such that its ground state wave function (and hence 
the boson masses fi a ) reproduces the same potential: 

V N (<p)^f(<p)^n a ->V N (<p). (4) 

3 Chiral dynamics 

It is straightforward to determine the chiral condensate. In the large or small 
volume limit (or equivalently in the low or high temperatue limit) an analytic 
expression can be obtained. For TV > 3 with degenerate fermion masses (m a = 
m <C m) 



1 ™ " \ «+t /m\ jv+t „ _ jv i 

-— 2e 7 cos— — for T < m«+i/i«+i 

4A m N) W (5) 

TT(N - 1) fl 

Here 9 is defined in the interval -tt < 9 < +ir by 9 = 9 - 2tt[(9 + tt)/2tt\. 

A few conclusions can be drawn here. First of all, at low T the chiral 
condensate is not analytic in m. This^point was noticed in the N = 2 case at 
T = by Coleman twenty years agoQ Secondly, at T = 0, a cusp singularity 
appears at 9 = tt. Thirdly at high temperature perturbation theory in the 
fermion masses is applicable since the condensate becomes analytic in m. 
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Figure 1: ^-dependence of the chiral condensate at various temperatures in the N = 3 model 
with m/fi = .01. At T = 0, a cusp singularity appears at 8 = 7r. 

At moderate temperatures Eq. (0) must be solved numerically. In fig. 1 
we have displayed the 9 dependence of chiral condensate in the N = 3 case 
with m/fx = 0.01 at various temperatures. One can see the cusp singularity 
develop as T approaches zero. 

At very low T (large L) the potential term dominates over the kinetic 
energy term in Eq. (g). In other words, the ground state wave function f(ip) 
has a sharp peak around the location of the absolute minimum of the potential 
Vn(<p)- As 9 varies, the location of the minimum also changes. In the case 
of degenerate fermion masses, the minimum is located at tp a — 9/N. The 
location of the minimum discontinuously shifts at 9 = 7r (mod 2n), which is 
the origin of the cusp singularity encountered in the 9 dependence of the chiral 
condensate at T = 0. 

When the fermion masses are not degenerate, the coefficients m a A a are 
all distinct. In the three flavor (N = 3) case the potential takes the form 

Va[tp] = — |gi cos^i + q 2 cos ^2 + q 3 cos(# — (fx — ^2)} (6) 

where all g Q 's are different. This potential has the same structure as the 
effective chiral Lagrangian in QCD. In the effective Lagrangian written by 



4 



WittenU a potential term reads 

V WiUcn {u) = ; 2| _ hft M( JJ + tft) + _A_(_ llndet U _ Of j (7) 

Here U is the pseudoscalar field matrix, whereas M = diag (m u ,md,m s ) is 
the quark mass matrix . The second term represents the contribution from 
instantons. The coefficient k is 0(1) in the large iV co i or limit. 

Phenomenologically m^, ^> m^m^, m^, which implies that k/N c ^> m a , 
or that upon diagonalizing U = diag (e^ 1 , e t4 > 2 , e^ 3 ), £>a = Hence (0) 
reduces to (^). These equations are not exactly the same however, as q a in (3) 
is not simply m a but depends on {m a } rather nontrivially. 

In the past various conclusions were drawn based on (Q), which turn out 
to be perfectly correct in our context as well. Physical quantities are periodic 
in 9 with period 27r. With degenerate quark masses a cusp singularity appears 
at 8 = n. Furthermore sufficiently large asymmetry in quark masses removes 
the singularity! 

Indeed, we have observed that at T — the location of the minimum of 
the potential determines the vacuum wave function. When q\ = q2 <C 93, 
i.e. the strange quark is heavy but the up and down quarks are degenerate 
(m u = trid -C m s ), we find that the discontinuous jump at 6 — ir in the location 
of the minimum remains. However, if one adds small asymmetry in the up and 
down quark masses (m u < rrid -C m s ), the minimum moves continuously to 
make a loop as 9 changes from — n to +n. In other words the cusp singularity 
disappears. 

4 Confinement 

One common method of examining confinement is to evaluate the Polyakov 
loop. In our method this is simply the expectation value of the Wilson line 
phase: 



Pn 



for — ^ integer 

(8) 

e"^ fMf(<Pa)*f(<P a + *&) for ! =n 

The result for q = e is displayed in fig. 2. Although P q shows a crossover 
transition and becomes vanishingly small at low T, the free energy F q of a test 
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Figure 2: T-dependence of the Polyakov loop and free energy in the N = 3 model with 
m/n = .01. The free energy stays finite at all T. 



charge q remains finite at all T . P q vanishes for a fractional charge. However, 
one cannot conclude that confinement of fractional charge from this result 
alone, as the vanishing of P q follows solely from gauge invariance. 

To get more information about confinement, one should insert a pair of 
sources, one with charge q and the other with — q, and examine the increase 
(or decrease) of the energy. The shift in the energy is parametrised as AE = 
ad + ■ ■ ■ where d is the distance between the two sources, and a is the string 
tension. 

In the multi-flavor case perturbation theory in mass cannot be employed. 
Nevertheless one arrives at a simple result. External charges are completely 
screened, but the effective 9 value is shifted between the two sources by an 
amount 2irq/e. This shift in turn changes. ±h.e chiral condensate and therefore 
the energy density between the sources&BEffl 

One finds 



Nm^(i/j t/j) g alt - (ipip)ej 



OeS = 1 



2irq 

e 



In particular, at T = 



N 
2~- 



n( 2e n) ( C0S Ar) -( cos iv) 



2N 
N + l 



2N ■ 
N + l 



(9) 



(10) 
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Notice that the string tension vanishes for an integer q/e. The string tension a 
is non-vanishing only when fcrmions are massive and chiral condensates have 
non-trivial 9 dependence. 

5 Heavy fermions 

When fermion masses become large, we need to solve a more general problem. 
For one flavor there is no (f degree of freedom. The vacuum wave function is 
expressed as f(pw) which must satisfy 



Here to = irfiL and k depends on the fermion mass and the wave function itself. 
One finds that 



The chiral condensate becomes bigger as the fermion mass m grows. The 
effect of a very heavy fermion never disappears even if m approaches infinity. 
Does this contradict the decoupling theorem, that heavy fermions are irrelevant 
in low energy physics? 

The answer is 'no'. Chiral condensates of heavy fields are not good mea- 
sures of low energy physics. Instead one should examine, for instance, the 
string tension defined when a pair of external sources is inserted as was done 
in the previous section. 

Eq. (JoJ) shows that the string tension a depends on two factors: to and the 

9 dependence of (ipip). Although m(ipip) becomes very large as to increases, 
its 9 dependence diminishes rapidly for a large to. In fig. 3 we have displayed 
the to dependence of a for 9 = and q/e — \. One can see that a increases 
for moderate to, but quickly approaches zero as m becomes large. Extremely 
heavy fermions are thus irrelevant for low energy physics. 

6 Summary 

There are similarities and differences between QCD4 and QED2: both have 
confinement, and their chiral dynamics are pretty much the same. However, in 

QED2 a non- vanishing string tension (ie. confinement) results only if m^i^) 




(11) 




— e 7 cos 9 for m -C jj, 



for to ^> /1 . 



(12) 
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Figure 3: Mass dependence of the string tension <r in the N = 3 model, a vanishes for a 
large fermion mass, which is consistent with the decoupling theorem. 



has non-trivial 9 dependence. In other words, if there were no chiral conden- 
sates, there would be no confinement in QED 2 . 

This can hardly be true in QCD4 where it seems that confinement and 
spontaneous chiral symmetry breaking are two separate phenomena. Chiral 
symmetry is spontaneously broken even in the chiral limit m a = 0, which is 
not the case in multiflavor QED2. Nonetheless lattice simulations, show that 
confinement and chiral symmetry breaking are intimately relatedcj 

Despite these differences, we can still learn quite a bit of QCD physics from 
QED2. After all we have determined various physical quantities such as chiral 
condensates, Polyakov loop, and string tension at any temperature and with 
arbitrary fermion masses. This is the luxury of two-dimensional gauge theory. 
In the work summarized here, we have presented a powerful method for solving 
QED2 which compliments lattice gauge theory and light-front methodstLil 
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